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ONE  DIMENSIONAL  D.C,  RESISTIVITY  DUE  TO 
STRONG  LOW  FREQUENCY  TURBULENCE 


The  problem  of  anomalous  plasma  d.c.  resistivity  in  one  dimension  is  one  of  the  oldest 
basic  plasma  physics  questions.  Its  resolution  has  eluded  the  plasma  community  despite  the  in- 
tensive theoretical  and  experimental  efforts''^.  Besides  its  fundamental  nature,  the  question  of 
anomalous  resistivity  is  of  utmost  practical  importance  for  laboratory  plasma  heating  and  for 
energy  dissipation  and  particle  acceleration  in  space  plasmas  According  to  our  present  under- 
standing, the  presence  of  a field  aligned  electron  current  in  the  plasma  with  a drift  velocity  w p, 
higher  than  a threshold,  excites  density  fluctuations  which  then  scatter  the  electrons  with  a 
collision  frequency  i/  * larger  than  the  small  angle  Coulomb  frequency  v,  thereby  producing  a 
higher  dissipation  rate.  This  provides  a consistent  and  experimentally  satisfactory  picture  if 
V*  > ilg  where  is  the  electron  cyclotron  frequency.  However,  when  i-*  < the  elec- 
tron magnetic  moment  is  conserved  and  the  electrons  should  be  treated  as  one  dimensional. 
In  this  case,  as  clearly  shown  by  Petviashvili,^  the  formation  of  a plateau  in  the  electron  distri- 
bution function  precludes  any  resistivity  enhancement  beyond  a few  percent  of  the  classical 
value.  It  is  the  purpose  of  this  letter  to  demonstrate  that  inclusion  of  strong  turbulence  effects 
can  remove  most  of  the  fundamental  difficulties  encountered  in  one-dimensional  anomalous 
resistivity  and  produce  a consistent  and  satisfactory  picture.  The  basic  idea  is  that  density  cavi- 
ties in  the  ion  background  cause  the  existence  of  low  frequency  potential  barriers.  When  a d.c. 
electric  field  is  applied  these  barriers  can  prevent  the  free  streaming  acceleration  of  a 

significant  fraction  of  the  electrons. 
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An  important  step  forward  in  understanding  when  anomalous  d.c.  resistivity  can  appear 
was  the  observation  of  Papadopoulos  and  Coffey'*  that  the  low  frequency  density  fluctuations 
could  be  generated  by  means  other  than  current  driven  instabilities.  This  removes  the  thres- 
hold requirement  on  v^  of  the  ambient  plasma.  In  their  paper  the  low  frequency  fluctuations 
appeared  due  to  the  presence  of  suprathermal  electron  beams.  Such  situations  exist  both  in  re- 
lativistic beam  plasma  heating  and  in  the  auroras  due  to  the  presence  of  the  energetic  precipi- 
tating electrons.  It  was  noted  that  the  presence  of  the  beams  creates  electron  plasma  oscilla- 
tions whose  pondermotive  force,  acting  on  the  plasma,  drives  low  frequency  density  fluctua- 
tions. In  going  into  the  one  dimensional  situation  and  within  the  framework  of  weak  turbu- 
lence theory  the  same  difliculties  reappear  as  in  the  current  driven  case.  However  under 

strong  turbulence  conditions  (i.e.,  — = > the  density  fluctuations  have  been  shown 

nT 

to  form  large  localized  density  cavities  (cavitons)  in  which  the  high  frequency  held  is  trapped 
(solitons)*'^.  Figure  1 shows  a typical  structure  as  seen  in  a computer  simulation  of  such  cavi- 
ties produced  by  the  high  frequency  (w^)  waves  due  to  the  existence  of  a beam. 

While  our  original  work  centered  on  such  beam-generated  cavities  which  occur  in  many 
situations  of  current  physical  interest,  the  basic  physical  effects  noted  in  this  letter  will  appear 
whenever  finite  amplitude  ion  cavities  exist  in  the  ambient  plasma.  A description  of  the  in- 
teraction of  the  background  electrons  with  the  cavities  requires  a local  theory  going  beyond  the 
limit  of  weak-turbulence  quasi-linear  wave-particle  interaction.  In  view  of  the  analytical  com- 
plexity of  the  subject  we  present  below  a phenomenological  theory  supplemented  by  computer 
simulations. 

THEORY 

Assume  that  there  is  a density  cavity  in  the  ion  background  with  a depth  ( n„  is 

the  average  ion  density.  The  local  density  is  n,  (.x)  — 
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rig).  In  order  to  maintain 
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charge  neutrality  at  low  frequency,  there  must  exist  a low  frequency  potential  internal  to  the 
plasma,  6<6,  large  enough  to  exclude  a fraction  of  the  electrons  equal  to  This  means  that 
inside  the  cavity 
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From  the  Vlasov  equation,  since /is  constant  along  lines  of  constant  energy  we  can  determine 


Jf,„  in  terms  of  If  is  a nondrifting  Maxwellian  then  (1  - ^njn^)  - e where 
8</»  is  normalized  to  the  electron  thermal  energy  and  v,^  = 1.  Defining  a trapping  velocity 
v,^,  = the  fraction  of  trapped  electrons  is 


This  is  plotted  on  Fig.  2 as  a function  of  cavity  depth.  In  the  absence  of  any  cavities,  the  elec- 


trons will  free  stream  with  an  average  velocity  Er  when  an  external  d.c.  field  is  ap- 

plied.  The  presence  of  the  density  cavities  divides  the  electrons  into  two  classes.  The  central 
part  of  the  distribution  with  |v|  < v„  does  not  accelerate,  since  it  is  reflected  by  the  potential 
barriers.  The  tail  with  initial  velocities  parallel  to  v^j  being  untrapped  is  displaced  by  v^^;  un- 
trapped particles  that  have  an  initial  velocity  opposite  to  v^^.  can  only  slow  down  to  a velocity 
''mm  = (8<<>)''^.  At  that  point  they  are  scattered  by  180°  by  the  potential.  This  means  that 
the  electrons  in  the  left-hand  tail  have  undergone  a sudden  increase  in  velocity  equal  to  2 
v^j„.  Therefore,  the  average  electron  drift  velocity  is 
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SIMULATION  RESULTS 


(2) 


We  present  here  a series  of  ID  computer  simulations  showing  the  acceleration  of  a plas- 
ma by  a constant  d.c.  electric  field.  A particle-fluid  hybrid  code’  and  a Vlasov  simulation 
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code^  were  used.  In  the  first  group  of  simulations  we  imposed  a fixed  ion  density  fluctuation 
on  the  plasma.  This  allowed  us  to  vary  easily  ^njng  and  the  shape  of  the  fluctuation  in  order 
to  assess  their  effect.  In  the  second  group  of  simulations,  the  ion  fluctuations  were  generated 
self  consistently  by  the  pondermotive  force  of  a high  frequency  long  wavelength  pump,  such 
as  expected  in  the  presence  of  weak  electron  beams  or  a laser. 

Figure  3 shows  the  results  of  the  first  set  of  runs.  The  straight  line  is  the  free  streaming 
velocity  of  the  plasma  and  occurs  when  6/7,  —0.  The  curved  lines  are  calculated  from  equation 
2 assuming  — v,,.  Three  of  the  simulations  plotted  used  the  particle  fluid  hybrid  code. 
The  electrons  were  particles.  The  ions  were  a fixed  fluid  background.  The  small  differences  of 
the  simulations  with  the  theoretical  curves  can  be  attributed  to  the  electron  noise  due  to  the 
finite  number  of  particles.  As  can  be  seen  in  Fig.  3 this  was  actually  confirmed  by  using  the 
noise  free  Vlasov  code  with  /w,  — 00.  This  code  due  to  its  noise  free  property  allowed  us  to 
clearly  distinguish  between  trapped  and  untrapped  electrons.  Figure  4 shows  the  electron  dis- 
tribution at  r - 0 and  - 750.  With  - .2,  50%  of  the  particles  were  trapped  in 

agreement  with  the  theory.  Whitfield  and  Skarsgard’  observed  this  splitting  of  the  electron 
distribution  function  in  particle  simulations  in  which  they  modelled  the  effect  of  a bumpy  mag- 
netic field  with  a fixed  sinusoidal  density  variation. 

We  initialized  the  ions  in  the  Vlasov  simulation  four  different  ways 
(n,  (x)  - flg  (1 -f  (8rt,/rt„)  sin  fcf) );  (1)  A:  — 2rr/L  (2)  A:  - ir/L  (3)  A:  - tr/L,  n,  (x)  — 
for  X < L/2.  (4)  k - w/L,  n.  - for  x < 3 £72  (only  a cavity).  L,  the  grid  length,  equals  50 
\£).  8n,/ng  — 0.2  in  all  four  cases.  It  was  found  that  the  interaction  was  independent  of  both 
the  width  and  the  spacing  of  the  cavities  and  dependent  only  on  the  value  of  the  density 
minimum.  These  results  are  consistent  with  the  simplified  theory  presented  previously. 
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We  examined  next  two  situations  where  the  ion  cavities  are  self  consistently  generated 
due  to  the  OTSl. 

Figure  5 shows  the  results  of  a hybrid  simulation  with  w,  — 128w,  where  a kinetic  beam 
plasma  instability  with  growth  rate  y — .006tu^,  and  phase  velocity  — 81v,^  was  driven. 
The  mode  was  nonlinearly  stabilized^  by  creating  a set  of  10  density  cavities  (see  Figure  1). 
The  rms  value  of  the  density  fluctuations  was  15%,  while  the  localized  cavity  depth  was  .2-.3. 
An  electric  field  was  turned-on  at  (o^r  — 700.  It  can  be  seen  that  the  value  of  the  accelera- 
tion was  substantially  smaller  than  free  streaming. 

In  another  run  the  Vlasov  code  was  used  with  a constant  dipole  pump  of  energy 
£2 

— — .1.  Such  a run  simulates  the  possibility  of  a microwave  or  laser  pump  at  a>  across 

the  current  carrying  system.  A large  single  cavity  was  formed  in  this  case.  The  plasma 
response  was  again  consistent  with  the  simple  theoretical  picture. 

SUMMARY  AND  CONCLUSIONS 

The  results  of  this  letter  should  be  considered  as  a preliminary  but  major  step  in  under- 
standing and  modeling  the  long  times  anomalous  resistivity  in  one  dimension.  Referring  first 
to  our  self  consistent  results  of  the  second  simulation  group  we  demonstrated  that 

(a)  In  the  presence  of  electron  beams  or  externally  imposed  high  frequency  (ai^)  elec- 
tromagnetic fields,  large  local  density  cavities  (Sn/n^  =;  .2  —.6)  can  be  formed. 

(b)  The  potentials  associated  with  these  cavities  are  large  enough  to  trap  a significant 
fraction  of  the  thermal  electrons,  so  that  the  entire  current  will  be  carried  by  a small  fraction 
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of  untrapped  electrons.  This  effect  will  experimentally  appear  as  an  enhanced  resistance  and 
the  appearance  of  localized  runaway  beams.  The  long  time  scale  behavior  of  the  system  will 
depend  on  how  these  beams  are  thermalized.  A possible  mechanism  is  further  two  stream  in- 
teraction of  the  runaways  with  the  trapped  electrons.  This  process  will  depend  critically  on  the 
spacing  among  cavities  and  is  presently  under  study. 

Furthermore  the  results  of  the  first  simulations  demonstrated  that  similar  resistivity  will 
appear  anytime  an  instability  can  create  finite  amplitude  density  fluctuations.  The  ion  cyclo- 
tron instability  seems  to  be  a good  candidate  and  is  currently  under  study.  In  the  presence  of 
an  electron  current,  the  generation  of  ion  cavities  will  disrupt  the  current.  By  recalculating  Eq. 

(2)  using  a drifting  Maxwellian  one  obtains  a formula  in  agreement  with  preliminary  computer 
simulations. 

In  concluding  we  mention  that  these  concepts  can  help  in  understanding  return  current 
heating  in  beam  plasma  interactions  and  extend  the  Papadopoulos-Coffey  anomalous  resistivity 
in  the  auroral  zones  to  regions  where  = 1.  The  results  of  the  constant  pump  runs 

suggest  the  possibility  that  a laser  created  corona  with  large  cavities  a pellet,  can  produce  short 
electron  beam  deposition  lengths  for  e-beam-pellet  fusion. 
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